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ABSTRACT 

Context. It appears that most (if not all) massive stars are born in multiple systems. At the same time, the most massive binaries are 
hard to find due to their low numbers throughout the Galaxy and the implied large distances and extinctions. 

Aims. We want to study: [a] LS III -1-46 11, identified in this paper as a very massive binary; [b] another nearby massive system, 
LS III +46 12; and [c] the surrounding stellar cluster, Berkeley 90. 

Methods. Most of the data used in this paper are multi-epoch high-S/N optical spectra though we also use Lucky Imaging and archival 
photometry. The spectra are reduced with devoted pipelines and processed with our own software, such as a spectroscopic-orbit code, 
CHORIZOS, and MGB. 

Results. LS III+46 11 is identified as a new very-early-O-type spectroscopic binary [03.5 If* + 03.5 If*] and LS III+46 12 
as another early 0-type system [04.5 V((f))]. We measure a 97.2-day period for LS III +46 11 and derive minimum masses of 
38.80+0.83 Mq and 35.60+0.77 Mq for its two stars. We measure the extinction to both stars, estimate the distance, search for optical 
companions, and study the surrounding cluster. In doing so, a variable extinction is found as well as discrepant results for the distance. 
We discuss possible explanations and suggest that LS III +46 12 may be a hidden binary system, where the companion is currently 
undetected. 

Key words. Binaries: spectroscopic — Dust, extinction — Open clusters and associations: individual: Berkeley 90 — Stars: early- 
type — Stars: individual: LS III +46 11 — Stars: individual: LS III +46 12 


1. Introduction 

Multiplicity is an endemic disease among O stars 
(iMason et al.lll998l:ISana et aT]l2013atlSota et al.ll201^ . Indeed, 
it i s difficult to find a n O star that was born in isolation; among 
the ISota et al.l (l2014l) sample of 07-09.7 V-IV stars south of - 
20 deg, there is just one clear, currently single object, namely 
p Col, which is however a known runaway and thus likely was 
born in a multiple system or compact cl uster. Only a smal l frac¬ 
tion of the other types of O stars in the ISota et al] (|2014|) sam¬ 
ple are apparently single but those cases are more luminous and 


can easily hide an e.g. main-sequence B-type companion in their 
glare. Both types of binarity, spectroscopic and visual, are com¬ 
mon among O stars and combinations of both types (implying 
higher-order multiplicities) are also relatively frequent, usually 
in the form of hierarchical systems, with one pair in a close or- 
bit and additional s tar(s) at larger separations. The OWN survey 
(iBarba et al.l 120101) has detected a peak around 10 days in the 
spectroscopic period distribution of 240 southern massive stars. 
Most of those systems will i nteract during the ir evolution and 
many mergers are expected (ISana et al.ll2012l) . Even in those 
cases where interactions will not take place it is important to 
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LS III +46 11 

LS III +46 12 

Reference 

Sp. type 

03.5 If* + 03.5 If* 

04.5 V((f)) 

This work 

RA (J2000) 

20:35:12.642 

20:35:18.566 

H0g et al. 12000) 

dec (J2000) 

+46:51:12.12 

+46:50:02.90 

Hpg et al. 12000) 

1 (deg.) 

84.8844 

84.8791 

H0e et al. 12000) 

b (deg.) 

+3.8086 

+3.7836 

H0g et al. 12000) 

Ly 

10.889+0.021 

10.268+0.009 

This work 

Ks 

6.971+0.023 

7.470+0.023 

Skrutskie et al. 12006). this work 


Table 1. LS III -h 46 11 and LS III -h 46 12 summary. 



Fig. 1. 2MASS KsHJ three-color RGB mosaic of Berkeley 90. 
The intensity level in each channel is logarithmic. 


characterize the binarity of O stars since ignoring it can lead to 
biases in the measured properties, such as derived masses and 
predicted colors of small stellar populations or young clusters. 

The issue of multiplicity among massive stars also affects 
the measurement of the stellar upper mass limit, since one has 
to make sure that a very luminous object is not in reality a com¬ 
bination of two or more objects, an issue that has produced a 
numbe r of false alarms in the past, see e.g. lMafz Apellaniz et al.l 
(l2007h . For a long time, all very massive stars were thought to 
be born as very early O-type objects but the most current infor¬ 
mation lean s towards the extreme ly massive ones being born as 
WNh stars (ICrowther et al.ll^lOh . However, statistics are very 
poor to consider the case closed. For example, only one 02 star 
with accurate blue-violet spec tral classification is known in the 
Milky Way, HD 93 129 AaAb (ISota et al.ll2Q14li . and it is a mul¬ 
tiple system located in the center of a dense stellar cluster. We 
clearly need to find more very massive stars to understand their 
properties and evolution. 

The sparse young open cluste r Berkeley 90 (Fig.[TJ is almost 
unstudied, but ISandulea'5 (Il974l) suggested that the early-type 
stars LS III h- 46 11 and LS III -1-46 12 were members of the clus¬ 
ter. Both o bjects are classified as OB^ in the Luminous Stars 
catalogue (iHardom et al.l [T9641) . suggesting hi gh luminosities . 
From an analysis of 2MASS data for the regioni lTadrossI (l2008h 
concluded that the cluster was young (<100 Ma old), located 


at 2.4 kpc, and reddened by E(B - V) - 1.15 mag. IWendkeil 
(Il97lh suggested that LS III h-46 11 and LS III -1-46 12 were the 
ionizing sources for the Hit region [GSSSI 215, which can be 
identified with Sh 2-115 (iHarten & Fellill980h . lMaver & Mac^ 
(Il973h derived a spectral type 06 and a distance of 2.3 kpc for 
LS III h-46 12, but were unable to provide a spectral type for 
LS III h-46 11, which is more reddened. LS III h-46 11 was iden¬ 
tified as the counterpa rt of the moderately hard ROSAT source 
RX J2035.2H-4651 by iMotch et al.l (1 19971) . during a search for 
possible high-mass X-ray binaries. They derived a spectral type 
03-5 III(f)e from intermediate-resolution spectroscopy, and con¬ 
cluded that LS III h-46 11 could produce the observed X-ray lu¬ 
minosity of ~ 5 X 10^^ ergs ' with “probably no need” for a 
compact companion but that “the extreme value of the X-ray lu¬ 
minosity of this star would certainly deserve further detailed in¬ 
vestigation” without mentioning an alternative (e.g. a colliding- 
wind binary). iMotch et ^ (Il997h also detected an X-ray source 
coincident with LS III h-46 12 with about one third of the X-ray 
flux ofLS III h-46 11. 


As part of our Galactic O -St ar Sp ectroscopic Survey 
(GOSSS, iMafz Apellaniz et~^ 1201 ll) . we observed 
LS III h-46 11 with the TWIN spectrograph of the 3.5 m 
telescope at Calar Alto in early November 2009. We immedi¬ 
ately confirmed the very early nature of the system but more 
surprising was the double-lined spectroscopic binary (SB2) 
nature detected in the He ii lines. Given the importance of the 
discovery, we observed the system on three consecutive nights 
(1h-2h-3 November) and we noticed only a small decrease in the 
velocity separation, thus excluding a short period. For the next 
two years we reobserved LS III h-46 11 on different occasions 
but we failed to detect clearly split He ii lines until September 
2011, when we observed it with the 4.2 m WHT at La Palma. 
The new detection gave us a series of possible periods for the 
system and led to new observations in the subsequent years. By 
2013 we had a preliminary orbit with a period close to 97 days 
and by 2014 we obtained the final orbit, with a large eccentricity 
that explained the measurement difficulties. 


In this paper we first describe the used data, both the pri¬ 
mary spectroscopic observations and the complementary high- 
resolution imaging and archival photometry. We then present 
our results: the spectral classification of LS III h-46 11 and 
LS III h-46 12, the spectroscopic orbit of LS III h-46 11, the ex¬ 
tinction and distance to both stars, a search for visual compan¬ 
ions, and a global analysis of Berkeley 90. We finally discuss our 
results, including their relevance for the study of the stellar upper 
mass limit, and present our conclusions. In a subsequent paper 
we will analyze the foreground ISM in front of Berkeley 90. 


2 


















































J. Mafz Apellaniz et al.: The little-studied cluster Berkeley 90 


date 

(YYYY-MM-DD) 

HJD-2 400000 

Telescope 

Wavelength range 
(A) 

2009-11-01 

55 137.317 

CAHA-3.5 m 

3950-5050 

2009-11-02 

55 138.313 

CAHA-3.5 m 

3950-5050 

2009-11-03 

55 139.320 

CAHA-3.5 m 

3950-5050 

2009-11-03 

55 139.503 

CAHA-3.5 m 

3950-5050 

2009-12-05 

55 171.380 

OSN-1.5 m 

5350-6750 

2009-12-07 

55 173.355 

OSN-1.5 m 

5350-6750 

2009-12-09 

55 175.317 

OSN-1.5 m 

5350-6750 

2009-12-10 

55 176.363 

OSN-1.5 m 

5350-6750 

2010-05-01 

55 318.674 

OSN-1.5 m 

5350-6750 

2010-06-17 

55 365.659 

OSN-1.5 m 

5350-6750 

2010-06-26 

55 374.669 

OSN-1.5 m 

5350-6750 

2011-06-13 

55 726.513 

WHT 

3890-5570 

2011-06-15 

55 728.545 

WHT 

3890-5570 

2011-06-21 

55 734.768 

HET 

3811-4709 + 4758-5735 

2011-06-26 

55 739.779 

HET 

5311-6275 + 6396-7325 

2011-09-09 

55 814.597 

WHT 

3890-5570 

2011-09-10 

55 815.443 

WHT 

3890-5570 

2011-09-11 

55 816.460 

WHT 

3890-5570 

2011-09-12 

55 817.381 

WHT 

3890-5570 

2011-09-12 

55 817.455 

NOT 

3690-7180 

2011-09-12 

55 817.568 

WHT 

3890-5570 

2011-09-13 

55 818.791 

HET 

5311-6275 + 6396-7325 

2011-09-15 

55 820.489 

NOT 

3690-7180 

2011-09-17 

55 822.365 

NOT 

3690-7180 

2011-09-28 

55 833.751 

HET 

3811-4709 + 4758-5735 

2011-10-02 

55 837.724 

HET 

5311-6275 + 6396-7325 

2011-10-03 

55 838.724 

HET 

5311-6275 + 6396-7325 

2011-10-06 

55 841.706 

HET 

5311-6275 + 6396-7325 

2011-10-07 

55 842.706 

HET 

5311-6275 + 6396-7325 

2011-10-10 

55 845.719 

HET 

5311-6275 + 6396-7325 

2011-10-31 

55 866.649 

HET 

5311-6275 + 6396-7325 

2011-11-15 

55 881.601 

HET 

5311-6275 + 6396-7325 

2011-11-27 

55 893.590 

HET 

5311-6275 + 6396-7325 

2012-04-11 

56029.953 

HET 

5311-6275 + 6396-7325 

2012-05-15 

56063.877 

HET 

5311-6275 + 6396-7325 

2012-05-28 

56076.821 

HET 

5311-6275 + 6396-7325 

2012-06-10 

56089.815 

HET 

5311-6275 + 6396-7325 

2012-07-06 

56115.596 

WHT 

3910-5590 

2012-07-08 

56117.657 

WHT 

3910-5590 

2012-08-27 

56 167.808 

HET 

3811-4709 + 4758-5735 

2012-09-07 

56178.537 

CAHA-2.2 m 

3925-9225 

2012-09-23 

56 194.424 

WHT 

3910-5590 

2012-09-24 

56 195.442 

WHT 

3910-5590 

2012-09-25 

56 196.480 

WHT 

3910-5590 

2012-10-02 

56203.713 

HET 

5311-6275 + 6396-7325 

2012-10-04 

56205.732 

HET 

5311-6275 + 6396-7325 

2012-10-05 

56 206.727 

HET 

5311-6275 + 6396-7325 

2012-10-22 

56223.343 

CAHA-2.2 m 

3925-9225 

2012-12-04 

56266.548 

HET 

5311-6275 + 6396-7325 

2012-12-05 

56267.377 

CAHA-2.2 m 

3925-9225 

2012-12-07 

56269.531 

HET 

5311-6275 + 6396-7325 

2012-12-30 

56292.280 

CAHA-2.2 m 

3925-9225 

2013-06-14 

56458.541 

WHT 

3910-5590 

2013-09-30 

56566.341 

CAHA-2.2 m 

3925-9225 

2013-10-15 

56581.411 

WHT 

3910-5590 

2013-10-16 

56 582.442 

WHT 

3910-5590 

2013-10-17 

56583.344 

WHT 

3910-5590 

2014-04-13 

56761.624 

CAHA-2.2 m 

3925-9225 

2014-04-14 

56762.611 

CAHA-2.2 m 

3925-9225 

2014-06-07 

56816.568 

CAHA-3.5 m 

3880-4980 

2014-07-12 

56851.562 

CAHA-3.5 m 

3950-5050 

2014-08-04 

56874.526 

NOT 

3690-7180 

2014-08-05 

56875.513 

CAHA-2.2 m 

3925-9225 

2014-08-05 

56875.572 

NOT 

3690-7180 

2014-08-06 

56876.373 

CAHA-2.2 m 

3925-9225 

2014-08-07 

56 877.446 

CAHA-2.2 m 

3925-9225 

2014-08-08 

56878.583 

NOT 

3690-7180 

2014-08-14 

56884.384 

CAHA-2.2 m 

3925-9225 

2014-08-16 

56886.528 

CAHA-2.2 m 

3925-9225 

2014-08-18 

56888.542 

CAHA-2.2 m 

3925-9225 


Table 2. Spectroscopic observation log for LS III -1-46 11. The 
date refers to the local evening. 


date 

(YYYY-MM-DD) 

HJD-2400000 

Telescope 

Wavelength range 
(A) 

2011-06-13 

55 726.513 

WHT 

3890-5570 

2011-06-15 

55 728.545 

WHT 

3890-5570 

2011-09-09 

55 814.597 

WHT 

3890-5570 

2011-09-10 

55 815.443 

WHT 

3890-5570 

2011-09-11 

55 816.460 

WHT 

3890-5570 

2011-09-12 

55 817.381 

WHT 

3890-5570 

2011-09-12 

55 817.568 

WHT 

3890-5570 

2012-07-06 

56115.596 

WHT 

3910-5590 

2012-07-08 

56117.657 

WHT 

3910-5590 

2012-08-15 

56 155.648 

HET 

5311-6275 + 6396-7325 

2012-09-23 

56 194.424 

WHT 

3910-5590 

2012-09-24 

56 195.442 

WHT 

3910-5590 

2012-09-25 

56 196.480 

WHT 

3910-5590 

2013-06-14 

56458.541 

WHT 

3910-5590 

2013-10-15 

56 581.411 

WHT 

3910-5590 

2013-10-16 

56 582.442 

WHT 

3910-5590 

2013-10-17 

56583.344 

WHT 

3910-5590 

2015-03-02 

57 084.715 

CAHA-2.2 m 

3925-9225 

2015-03-03 

57 085.717 

CAHA-2.2 m 

3925-9225 

2015-03-28 

57 110.696 

CAHA-2.2 m 

3925-9225 


Table 3. Spectroscopic observation log for LS III -1-46 12. The 
date refers to the local evening. 


2. Data 

2.1. Spectroscopy 

The spectroscopic data for LS III -1-46 11 used in this pa¬ 
per were obtained with six different telescopes and are listed in 
Table |3 Here we describe the data grouped under the different 
projects in which they were obtained. 

1. GOS SS: the Galactic O-St ar Spec troscopic Survey is de¬ 
scribed in iMalz Apellaniz et’aD (1201 Ih . The GOSSS spectra 
used here were obtained with the TWIN spectrograph of the 
3.5 m telescope at Calar Alto (CAHA) and the ISIS spectro¬ 
graph of the 4.2 m William Herschel Telescope (WHT) at La 
Palma. The spectral resolving power was measured in each cas43 
and determined to be between 2800 and 3400 (TWIN) and be¬ 
tween 3900 and 5400 (ISIS). The typical S/N per resolution el¬ 
ement of the data is 300. The spectral range of the WHT data 
includes both Heii 44541.591 and Heii 45411.53 while that of 
the CAHA-3.5 m data included only the formefl GOSSS is 
also obtaining spectroscopy with the Albireo spectrograph of 
the 1.5 m telescope at the Observatorio de Sierra Nevada (OSN) 
but its aperture is too small to obtain a good S/N in the blue- 
violet. However, we also obtained some OSN yellow-red spectra 
at ~ 2000 in order to study He ii 45411.53. 

2. N oMaDS: the Northern Mas sive Dim Survey is de - 
scribed in iMafz Apellaniz et al.l (1201 2h and lPellerin et alJ (1201 2h . 
NoMaDs spectra were obtained with the High Resolution 
Spectrograph of the 9-m Hobby-Eberly Telescope (HET) at 
McDonald Observatory. Their spectral resolving power is 30 000 
and all epochs include Heii 45411.53 though two different se¬ 
tups were used on different occasions: one that goes from the 
violet to the yellow with a small gap (3811-4709 A -i- 4758- 


’ Note that the GOSSS spectra shown in ISota et al.l ( 1201 ll . [20l^ have 
been homogenized to R = 2500 whenever the resolving power was 
slightly higher. 

^ Throughout this work we refer to each line with the exact rest wave¬ 
length we are using in order to avoid problems with future comparisons 
of absolute velocity values. The different number of significant digits 
reflects the uncertainty in our knowledge of the rest wavelengths and/or 
the single/multiplet nature 
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5735 A) and one that goes from the green to the red with another 
small gap (5311-6275 A -H 6396-7325 A). 

3. lACOB: the Instituto de Astrofisica d e Can arias 
OB database is describe d in ISimon-Dfaz et alJ (1201 111 and 
ISimon-Diaz et alJ (l2015bl) . lACOB spectra were obtained with 
the FIES spectrograph at the 2.6 m Nordic Optical Telescope 
(NOT) at la Palma. Their spectral resolving power is 23 000 and 
they cover the whole optical range. Note that some epochs were 
obtained as part of lACOB itself and some were obtained during 
service time after a request to observe LS III -1-46 11 was granted 
additional time. 

4. CAFE-BEANS: the Calar Alto Fiber-fed Echelle 
Binary Evolution Andalusian Northern Survey is described in 
iNegueruela et alJ (l2014l) . CAFE-BEANS spectra were obtained 
with the CAFE spectrograph at the 2.2 m telescope at Calar 
Alto. Their spectral resolving power is 65 000 and they cover 
the whole optical range. 

The spectroscopic data for LS III -1-46 12 were obtained with 
GOSSS, NoMaDS, and CAFE-BEANS and are listed in Table[3] 
Note that the GOSSS instruments are long-slit spectrographs so 
the LS III -1-46 11 and LS III -1-46 12 spectra were obtained si¬ 
multaneously. 

The G OSSS spectroscopy was reduced with the pipeline de¬ 
scribed by ISota & Mafz Anellanizl (1201 Ih while the NoMaDS, 
lACOB, and CAFE-BEANS were reduced with specific 
pipelines for each instrument. The telluric lines in the high- 
re solution_specQ'a were eliminated according to the procedure 
of iGardini et alJ (l2013h . Note that the high-resolution spectra 
cover a large range in wavelength and that LS III +46 11 and 
LS III -1-46 12 are moderately reddened. That means that the S/N 
is highly variable as a function of wavelength and that in the vi¬ 
olet extreme only the NoMaDS spectra have a good S/N. 


2.2. High-resolution imaging 


To complement the spectroscopic observations of 
LS III -1-46 11 and LS III -1-46 12 we also present high-resolution 
imaging data that are part of a visual multiplicity survey of 
massive stars we are conducting using the Lucky Imaging 
instruments AstraLux Norte (at the 2.2 m CAHA telescope) 
and AstraLux Sur (at the 3.6 m NTT telescope at La Silla). 
The Lucky Imaging technique involves taking a large number 
(typically 10000) of very short exposures (typically, 30 ms), 
selecting a small fraction of those based on the PSF quality, and 
then combining them with a drizzle-type algo rithm. The survey 
itself was presented in iMafz Anellanizl (l2010l) : the reader is re¬ 
ferred to that paper for details on the data and proces sing. Later 
AstraL ux re sults for individua l stars a re discussed in Sota et aki 
(l2014ll and ISimon-Dfaz et al.l (l2015ah . iMafz Anellanizl (l2010h 
used a Sloan-z-band image of LS III -1-46 11. Here we use that 
image plus two additional z-band AstraLux Norte images of 
LS III -1-46 11, two Sloan-/-band images of LS III -1-46 11, and 
one z-band image of LS III -1-46 12, obtained between 2007 and 
2013. 


2.3. Archival photometry 

We used Simbad and Vizier to compile photometric data 
(Johnson UBV, Tycho-2 BV, and 2MASS JHKs) from the liter¬ 
ature for LS III -1-46 11 and LS III -1-46 12. The photometry was 
processed with CHORIZOS, as explained in the next section. 

We also compiled the IPHAS and 2MASS photometry of 
the Berkeley 90 and foreground populations. For the case of 


LS III +46 11, the 2MASS PSC gives J and El magnitudes but 
not a value, where an X flag appears, indicating that “’there is 
a detection at this location, but no valid brightness estimate can 
be extracted using any algorithm”. For LS III +46 12, which has 
similar NIR magnitudes (see Table |8]l, the 2MASS PSC gives 
the three magnitudes with good quality flags. We downloaded 
the 2MASS images to look into this issue. 

LS III -1-46 11 and LS III -1-46 12 have NIR magnitudes that 
make them slightly saturated (by ~1 magnitude) in the 2MASS 
images (which have saturation levels of 9.0, 8.5, and 8.0 for 
J, H, and Ks, respectively). For that reason, their 2MASS PSC 
magnitudes were obtained by aperture photometry in the 51 ms 
“Read_l” exposures which, unfortunately, are not publicly avail¬ 
able. However, for such a small level of saturation only 1-4 pix¬ 
els (depending on the star centering) are expected to be satu¬ 
rated so it should be possible to do differential (between the two 
stars) aperture photometry with an inner and an outer integra¬ 
tion radius. We tested this on LS III +46 11 and LS III -1-46 12 
using an inner radius of 1.5 pixels and an outer radius of 4.0 pix¬ 
els and found A7 = 0.304 and A// = 0.412, values that are 
within one sigma of those in the 2MASS PSC (0.281+0.031 
and 0.424+0.025, respectively, see Table |8]l, lending credibil¬ 
ity to the technique. We looked at the LS III +46 11 Xs radial 
profile and noticed no peculiarities so we applied the same tech¬ 
nique, which yielded AXs = 0.499. Therefore, we derive a value 
of Xs = 6.971 + 0.023 for LS III +46 11, where we adopted as 
uncertainty the same value as for LS III +46 12. 


3. Results 

3.1. Spectral classification ofLS III +46 11 and LS III +46 12 

The top two plots of Figures|2]and[3show two GOSSS spec¬ 
tra of LS III +46 11, one near maximum velocity separation and 
one close to the point where both stars have the same velocity. 
In the first plot, the spectra of the two components are virtually 
indistinguishable, with all the relevant stellar lines showing very 
similar ratios and widths. The only difference can be seen in the 
different absolute depths, the ones of the primary being ~ 10% 
more intense than those of the secondary. Therefore, we assign 
the same spectral type to both components^ _ 

W e classified the spectra using MGB (iMafz Apellaniz et al.l 
l2Q12h and v2.0 of the GOSSS standard grid 
(iMafz Apellaniz et al.l 12014^ . In the middle plot of Figure |2] 
He I 44471.507 appears as a very weak line in comparison with 
Hen 44541.59 and Niv 44057.75 is in emission with a similar 
intensity as Niii 44640.64+46 41.85, indicating a spectral type 
of 03.5 (IWalborn et alJl200l . Hen 44685.71 has a P-Cygni 
prohl^l (which appears to be present in both components, see 
the top plot) yielding a luminosity class of I since II is not 
d ehned at 03.5. T herefore, given the required f suffix (Table 2 
in ISotaet al.ll20l4 . LS III +46 11 is an 03.5 If* + 03.5 If*. 

The bottom plot of Figure |2] shows a GOSSS spectrum of 
LS III +46 12. All the GOSSS spectra are consistent with a con¬ 
stant spectral type i.e. we detect no sign of the target being an 
SB2. The GOSSS spectra are not calibrated in absolute velocity 
(the spectra are left in the star reference frame) but a compari¬ 
son with prominent ISM lines shows no sign of relative velocity 
shifts, so we do not detect an SBl character for LS III +46 12, 
either. 

N IV 44057.75 is not seen in emission in the spec¬ 
trum of LS III+46 12 and He i 44471.507 is significantly 

^ Note that the absorption component is only slightly blueshifted. 
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Fig. 2. Sample GOSSS spectra of LS III -1-46 11 and LS III -1-46 12 in the classical blue-violet spectral classification range. The top 
spectrum shows an example of LS III h- 46 11 near maximum velocity separation, the middle one an example of LS III -1-46 11 near 
minimum velocity separation, and the bottom one an example of LS III h- 46 12. The three cases show WHT data with the original 
spectral resolving power. 


P 

To 

Kn 

e 

a> 

an sin i 

{Ml + M 2 ) sin^ i 

^epochs 

^dof 

tL 


97.193 ±0.010 days 
56002.504 ±0.088 HJD-2400000 
237.3 ± 1.2 km/s 

0.5685 ± 0.0036 
124.89 ± 0.62 degrees 

1.7438 ± 0.0077 AU 
74.9 ±1.0 Mo 

70 
65 
1.43 


Table 4. LS III ±46 11 spectroscopic binary results for Av s V 2 - 
vi. K \2 is the radial velocity amplitude of the relative orbit. 


stronger than in LS III ±46 11, yielding a spectral type of 
04.5. Hen 44685.71 is quite deep, so the luminosity class 
is V. Since Niii 44640.64±464L85 is clearly in emission. 
Cm 44647.419±4650.246±4651.473 is weak, and there are no 
signs of line broadening, the final spectral type is 04.5 V((f)). 



Fig. 3. Same as Fig. |2 for the wavelength range around 
He II 45411.53, the primary line we use for the stellar velocity 
measurements. 


Note that Niv 45200.41±5204.28±5205.15 is seen in ab- 
sorption in all cases, a s expected for these spectral types 
dOamen & Niemeill2002h . 
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p 

97.168 ±0.025 

days 

To 

56002.80 ±0.25 

HJD-2400000 


112.7 ±1.2 

km/s 

K 2 

122.9 ± 1.3 

km/s 

e 

0.5627 ± 0.0061 


Cl) 

126.1 ±1.2 

degrees 

71 

-17.68 ±0.98 

km/s 

72 

-20.88 ± 1.04 

km/s 

ai sin i 

0.8325 ± 0.0078 

AU 

a 2 sin i 

0.9073 + 0.0081 

AU 

M\ sin^ i7(Mi + Mj)^ 

8.15 ±0.98 

Mo 

Ml sin^ i/(Mi + M2f 

10.55 ± 0.28 

Mo 

M 2 IM 1 

0.9175 + 0.0095 


Ml sin^ i 

38.80 ±0.83 

Mo 

M 2 sin^ i 

35.60 ± 0.77 

Mo 

^epochs 

44 


^dof 

80 


tL 

0.95 



Table 5. LS III -i-46 11 spectroscopic binary results for vi + V 2 . 
Ki and K 2 are the radial velocity amplitudes of each orbit. 


3.2. The spectroscopic orbit of LS III +46 11 

Most spectroscopic orbits for OB stars are studied with He i 
lines (or even with metallic lines for B stars) because they are 
intrinsically narrower than Hen lines. However, for a system 
composed of two 03.5 stars He i lines are not practical because 
they are too weak and we are forced to resort to the broader 
He II line^. The strongest Hen optical absorption lines in an 
03.5 If* are Hen T4541.591 and Hen 45411.53. The large ex¬ 
tinction experienced by LS III -1-46 11 (see below) makes a given 
S/N easier to attain with the latter than with the former, so we se¬ 
lected He II /15411.53 as our primary line. Note, however, that the 
CAHA-3.5 m spectra do not include it. In those cases we used 
He II 44541.591 instead to measure the Av = V 2 -v\, the velocity 
difference between the primary and the secondary. Also note that 
the CAHA-3.5 m and OSN-1.5 m spectra were not calibrated in 
absolute velocity and that the WHT spectra were calibrated us¬ 
ing ISM lines (instead of the lamps) present in both the WHT and 
the high-resolution spectra (which were calibrated using lamps). 

All the epochs were htted simultaneously with an IDL code 
by leaving the intensity and width of each component fixed (al¬ 
lowing for the different resolutions of each spectrograph) and fit¬ 
ting a four-component Gaussian for He ii 45411.53 (the two stel¬ 
lar components plus the two DIBs at 5404.56 A and 5418.87 A, 
with the DIBs fixed at the same velocity for all epochs) and a 
two-component Gaussian for He ii 44541.591. The initial period 
search was conducted using an IDL implementation of the in¬ 
formation entropy algorithm of ICincotta et al.l (Il995h . The orbit 
fitting itself (including the final period calculation) was done in¬ 
dependently by two of us: [a] J.M. A. using a code developed by 
himself with the help of a previous routine w ritten by R.C.G. and 
[b1 R.H .B. using an improved version of the lBertiau & GrobbenI 
(1 19691) code. Results were compared and found to be compati¬ 
ble, so only the hrst ones will be reported here. 

A first analysis was carried out with Av s V 2 - vi and all 
the observations in Table |2] Results are displayed in Table |4] 
and plotted in the left panel of Figure |4] In some cases with Av 
close to zero it is not possible to distinguish which component 
was the one with the larger velocity, since the two components 
are blended into a single Gaussian and, given the similar fluxes. 

This is a factor that likely contributed to the previous failure of the 
detection of the SB2 character of LS III -1-46 11. 


switching them around does not appreciably change the reduced 
or x^g^ of the fit. For those epochs with identification confu¬ 
sion we used a Av of zero with large error bars. The orbit ob¬ 
tained with Av has a good fit (x^gd~ a period consistent 

with the initial estimate, a relatively large eccentricity (as ex¬ 
pected), and a very large minimum system mass (uncorrected 
for inclination) of 74.9 Mq. 

A second analysis was carried out with the separate vi and V 2 
measurements. For this second analysis we excluded the CAHA- 
3.5 m spectra (for which we only had Hen 44541.591 and were 
not calibrated in absolute velocity), the OSN-L5 m spectra (also 
uncalibrated in absolute velocity and with identification con¬ 
fusion in all cases, since we always caught LS III -1-46 11 near 
Av = 0 with that configuration), and the WHT spectra with iden¬ 
tification confusion. Results are displayed in Table|5]and plotted 
in the right panel of Figure |4] The value of x^^^ is even better 
than for the first analysis and the results are overall compatible. 
The uncertainties of P, Tq, e, and to are slightly larger due to the 
exclusion of some epochs. The mass ratio is close to the inten¬ 
sity ratio of the He ii 45411.53 or He ii 44541.59 lines and in the 
same sense: the component with the stronger lines is the most 
massive. 

The values of yi and 72 in Table |5] indicate that the 
LS III -1-46 11 center-of-mass radial velocity is in the range be¬ 
tween -21 km/s and -17 km/s using the Hen 45411.53 line. 
Fitting Gaussians to the LS III -1-46 11 high-resolution spectra 
with better S/N and lowest velocity separation also yield values 
in that range. However, doing the same with other lines gives 
different values: -27+2 km/s for Hei 45875.65, -23+4 km/s for 
O III 45592.25, and -11+3 km/s for Hen 48236.79. These dis¬ 
crepancies are not uncommon when analyzing early-type O stars 
due to the intrinsic width of their lines and the effect of winds. 
We can conclude that the true center-of-mass radial velocity of 
LS III +46 11 is between -25 km/s and -15 km/s without being 
able to provide a more precise measurement at this time. 

We also analyzed the velocity of LS III +46 12 using the 
four high-resolution epochs listed in Table [3 We were unable 
to detect clear radial velocity variations at a level of 10 km/s 
or higher but with such a small number of epochs it is not pos¬ 
sible to ascertain the spectroscopic binarity of the system. We 
measured the radial velocity of the system using the same four 
spectral lines as for LS III +46 11 and obtained a smaller range 
between -16 km/s and -10 km/s. The better agreement between 
lines is possibly caused by the weaker winds in the dwarf com¬ 
pared to the supergiant pair. These results are consistent with 
LS III +46 12 being a spectroscopic single with the same center- 
of-mass radial velocity as LS III +46 11 as the stronger winds of 
the latter likely bias its measured velocity towards more negative 
(blueshifted) values. Also, the experience with LS III +46 11 
and other systems indicates that we should not exclude the pos¬ 
sibility of LS III +46 12 being a spectroscopic system with a pe¬ 
riod of months or longer. 

3.3. CHORIZOS analysis of LS III +46 11 and LS III +46 12 

We processed the Johnson+Tycho+2MASS photometry of 
LS III +46 11 and LS III +46 12 using the latest version of 
the CHORIZOS code (iMafz At)ellanizll2004t) to determine the 
amount and type of extinction and the distance to the stars. The 
CHORIZOS runs were executed with the following conditions: 

- We used the Milky Way grid of iMafz Anellanizl (12013bh . 

in which the two grid parameters are effective tempera¬ 
ture (Teff) and photometric luminosity class (LC). The latter 
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Fig. 4. Phased radial velocity curves for Av (left) and v\ + V 2 (right). The color code identifies the telescope. Note that the left plot 
includes all data points while the right plot excludes the two telescopes without accurate absolute velocity calibration (CAHA-3.5 
m and OSN-1.5 m). 


X (Angstroms) 

5000 10000 20000 



X (Angstroms) 

5000 10000 20000 



1/X (pm ’) 


Fig. 5. Best SED CHORJZOS fits for LS III -1-46 11 (left, luminosity class of 1.0 used) and LS III -1-46 12 (right, luminosity class of 
5.0 used). Blue data points are used for the input photometry (vertical error bars indicate photometric uncertainties, horizontal ones 
approximate filter extent) and green stars for the synthetic photometry. 


quantity is defined in an analogous way to the spectroscopic 
equivalent but instead of being discrete it is a continuous 
variable that varies from 0.0 (highest luminosity for that Tetr) 
to 5.5 (lowest luminosity for that Tss). Note that the range is 
selected in order to make objects with spectroscopic lumi¬ 
nosity class V (dwarfs) have LC=s5 and objects with spectral 


luminosity class I have LC* 1. For O stars the spectral ener gy 
distributions (SEDs) are TLUSTY (iLanz & Hubenvl2003h. 
The ex tinction laws were those of iMafz Anellaniz et al.l 
(l2014bl) . which are a single-family parameter with the type 
of extinction defined by 7 ^ 5495 . The amou nt of extinction 
is para meterized by £(4405 - 5495). See iMafz Anellaniz! 
(l2()13ah for their relationship with Ry and E{B-V) and why 
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Quantity 

LS III +46 11 

LS III +46 12 

reff(K) 

41300 

41900 

luminosity class 

2 

^red 

1.0 

5.0 

1.36 

1.42 

£5495 

3.303±0.058 

3.377+0.040 

£(4405 - 5495) (mag) 

1.653±0.020 

1.255+0.011 

Av (mag) 

5.475±0.037 

4.272+0.021 

Vyo (mag) 

5.414±0.021 

5.995+0.017 

log d (pc, min) 

3.308 

3.103 

log d (pc, used) 

3.407 

3.172 

log d (pc, max) 

3.525 

3.276 


Table 6 . Results of the CHORIZOS fits for LS III -h 46 11 and 
LS III -h46 12. 


those quantities are not good choices to characterize extinc¬ 
tion. 

Th e Teff-spectral type c onversion used is an adapted version 
of iMartins et al. ( 2005h that in cludes the spectral subtypes 
used by Sota et al. ( 20111 l2014t) . 

Teff and LC were fixed while R 5495 , £(4405 - 5495), and 
log d were left as free parameters. The values of the Teft were 
established from the used TefF-spectral type conversion (see 
previous point). For LC we explored a range of possible val¬ 
ues by doing CHORIZOS runs with 0.0, 0.5, 1.0, 1.5, and 
2.0 for LS III -1-46 11 and runs with 4.0, 4.5, 5.0, and 5.5 for 
LS III h-46 12. 


The CHORIZOS results are shown in Table and the best 
SEDs are plotted in Figure |5] along with the input and synthetic 
photometry. 


As expected, the different runs for a given star with differ¬ 
ent values of LC give nearly identical results £ 5495 , 

and £(4405 - 5495) and only differ in logii. This happens 
because the optical-i-NIR colors of O stars are nearly inde¬ 
pendent of luminositjO for a fixed Teft- Therefore, here we 
will concentrate on the results for LC=L0 (LS III h-46 11) 
and LC=5.0 (LS III h-46 12) and consider the other runs only 
when discussing the distance. 

The values of indicate that the fit is good, even for 
two cases such as these where the extinction is consider- 
abl e. We also ran alterna tive CHORIZOS execu tions using 
the ICardelli et al.l (1198^ and iFitzpatric^ (Il999h extinction 


laws. For ICardelli et al.l (Il989l) the were similar, as ex¬ 


pected for stars with moderate extinction with broad-band 
photometry and £ ^495 values close to the canonical 3.1. For 
iFitzpatri^ (1 1 999h thex^^d were significantly worse (by a fac¬ 
tor of »2). T herefore, these results are ano ther sign of the 
validity of the iMafz Apellaniz et"^ (l2014bh extinction laws 
for Galactic targets. 

The two stars show values of £5495 which are compatible be¬ 
tween them and only slightly larger than the canonical value 
of 3.1. Therefore, the same type of dust appears to lie be¬ 
tween each star and us and the grain size is typical for the 
Milky Way. 

The LS III h-46 11 extinction is significantly higher than that 
of LS III h-46 12, which explains the similar magnitudes of 
the two objects despite LS III h-46 11 being expected to be 
intrinsically more luminous and located at the same dis- 


^ The largest difference takes place in the band due to the wind 
contribution but that effect is not taken into account in the TLUSTY 
SEDs and is only expected to he of the order of 0.01 magnitudes for the 
cases of interest here. 


tance. This result also implies that there is significant dif¬ 
ferential extinction within the Berkeley 90 field. The value 
of £(4405 - 5495) for LS III -h 46 1 2 is close to th e E(B - 
V) = 1.15 result for Berkeley 90 of iTadrossI (l2008h . 

- The values for log d listed are the uncorrected CHORIZOS 
output and are equivalent to spectroscopic parallaxes. They 
do not take into account the fact that LS III h-46 11 is an 
SB2 system with two components with similar luminosi¬ 
ties while LS III h-46 12 is apparently single. Therefore, 
the log d values for LS III h-46 11 have to be increased by 
log V 2 = 0.151. This implies that the derived dis¬ 
tances for LS III h-46 11 and LS III h-46 12 are incompati¬ 
ble. Another way to look at the discrepancy is that if we 
use the extinction-corrected apparent Johnson V magnitude 
(Vjfl) for LS III -h46 11 and apply a 2.5 log2 = 0.753 
correction, we end up with two stars with values be¬ 
tween 6.1 and 6 . 2 , not brighter but actually slightly fainter 
than LS III h-46 12. In other words, if the three stars are at 
the same distance, we would require the two objects with 
spectroscopic class I to be fainter than the object with spec¬ 
troscopic class V. We analyze the distance issue later on. 


3.4. Visual multiplicity 

Table I 2 ] gives the separation, orientation and magnitude dif¬ 
ference between the brightest star in the AstraLux Norte im¬ 
ages (either LS III h-46 11 or LS III h-46 12) and the rest of the 
detected point sources. Note that the AstraLux images are not 
absolutely calibrated, so only differential photometry is pro¬ 
vided; hence, all the data refer to stellar pairs. Nevertheless, 
the photometry in the i and z bands for the two bright stars 
can be derived from the best SED in the previous subsection: 
LS III h-46 11 a has zero-point-corrected AB i and z magni¬ 
tudes of 9.549 and 9.017, respectively, and LS III h-46 12 A has 
zero-point-corrected AB i and z magnitudes of 9.371 and 9.043, 
respectively. Table [ 8 ] gives the 2MASS JHKs magnitudes for 
the stars in the AstraLux images. All are detected except for 
LS III h-46 11 F, the dimmest star in the AstraLux image. Note 
that the 2MASS magnitudes are not complete, as usual for mod¬ 
erately crowded fields such as this one. 

The most relevant result is that LS III h-46 11 has no visual 
companions within 11" and that LS III h-46 12 has just one dim 
companion 6 " away. From the photometric point of view, this 
means that the analysis in the previous subsection does not ap¬ 
pear to include additional stars, so one can refer to the pho¬ 
tometry of LS III h-46 11 and LS III h-46 11 A as indistinguish- 
abl 43 (as we have done in the previous paragraph) and the same 
can be said about LS III h-46 12 and LS III h-46 12 A. From the 
physical point of view, this means that (barring any undetected 
components) LS III h-46 11 is a double system but likely not 
a higher-order one, since F or C (the closest companions) are 
likely to be unbound, especially considering that their environ¬ 
ment is the center of a cluster. On the other hand, LS III h-46 12 
B is ~ 10 000 AU in the plane of the sky awa y from A and 
could possibly be bound (iMalz Apellaniz! I 2 OIOI) . There is also 
another component visible in the 2MASS images (but not listed 
in Tables|2]or[8]because it fell just outside the 25"x25" AstraLux 


® LS III h-46 11 a includes two spectroscopic components but their 
separation, as previously derived, is of the order of 1 mas, two orders of 
m agnitude below what ca n be resolved with AstraLux Norte. See Fig. 2 
in iMafz Anellanizl i201(f) for the values of separations and magnitude 
differences that AstraLux Norte can resolve. 
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Pair 

Separation 

(") 

Orientation 

(degrees) 

Ai 

(mag) 

Az 

(mag) 

LS III +46 11 AB 

22.601+0.054 

139.94+ 0.18 

5.263+0.011 

5.209+0.012 

LS III +46 11 AC 

14.274+0.061 

114.23+ 0.15 

6.442+0.017 

6.260+0.014 

LS III +46 11 AD 

23.803+0.023 

131.77+ 0.07 

6.888+0.037 

6.445+0.040 

LS III +46 11 AE 

17.446+0.039 

123.09+ 0.34 

8.123+0.062 

7.571+0.036 

LS III +46 11 AF 

12.300+0.056 

105.54+ 0.22 

8.381+0.266 

7.802+0.226 

LS III +46 12 AB 

6.114+0.050 

318.83+ 0.20 

— 

6.968+0.034 

LS III +46 12 AC 

11.135+0.050 

318.60+ 0.20 

— 

7.781+0.053 


Table 7. AstraLux results for LS III +46 11 and LS III -1-46 12. 


Star 

2MASS ID 

J 

(mag) 

H 

(mag) 

K, 

(mag) 

Flag 

LS III +46 11 A 

J2035I264+465II2I 

7.653+0.023 

7.194+0.018 

6.971+0.023 

AAT 

LS III +46 11 B 

J20351402+4650549 

12.326 

12.050+0.054 

11.821+0.053 

UAA 

LS III +46 11 C 

J20351389+4651065 

13.402+0.049 

12.847+0.066 

12.585+0.070 

AAE 

LS III +46 11 D 

J20351436+4650562 

12.614 

12.638+0.033 

12.235+0.030 

UAA 

LS III +46 II E 
LS III +46 11 F 

J20351405+4651025 
undetected 

14.295+0.062 

13.373+0.049 

12.990+0.050 

AAA 

LS III +46 12 A 

J20351857+4650028 

7.934+0.021 

7.618+0.017 

7.470+0.023 

AAA 

LS III +46 12 B 

J20351823+4650072 

10.333 

10.158 

12.551+0.205 

UUC 

LS III +46 12 C 

J20351786+4650112 

14.450+0.124 

13.466+0.134 

13.242+0.050 

BBA 


Table 8. 2MASS photometry for the targets detected in the AstraLux images. ABC flags are used for normal detections (with 
increasingly larger uncertainties), U flags for upper limits on magnitudes, and T for measurements done in this work. 


Norte field of view) 7" to the E of LS III -1-46 12 A, thus increas¬ 
ing the probability of the existence of a bound companion. 

3.5. Berkeley 90 photometry 

The ar ea of Berkeley 9 0 is immersed in bright nebulosity 
(Sh 2-115. 1^ arr)lesslll959l) . indicating the existence of sources 
with large ionizing fluxes, but the cluster has surprisingly re¬ 
ceived ve ry little attention. We h ave used photo metric data from 
2MASS (ISkrutskie et al.l l2006h and IPHAS (iBarentsen et al.l 
l2014ll to study its properties. We have selected 2MASS sources 
within 3' of the nominal center of the cluster, as given in 
SIMBAD. We have rejected stars with bad quality flags, and 
cross-matched our selection with the IPHAS catalogue. As the 
IPHAS catalogue contains many spurious sources in areas of 
bright nebulosity, we have only accepted sources with a 2MASS 
counterpart within a (X'h radius. The Qir parameter, defined as 
Gir = {J - H) - 1.^ X (H - Ks), is very effective at separat¬ 
ing early and late-type stars (e.g.. IComeron & Pasaualil l2005t 
iNegueruela & Schurchl200^ . with early-type stars showing val¬ 
ues ai 0.0. We select objects with Qm < 0.08 as candidate early- 
type stars (see INegueruela & Schurchll200'^ . though emission¬ 
line stars also display negative values due to their excesses. 
The candidates are clearly concentrated towards the cluster cen¬ 
ter, confirming that they mainly represent the cluster population. 
The Ksl(J - Kg,) diagram for the resulting selection is displayed 
in Fig. 0 

Possible cluster members show a broad distribution in (J - 
Kg). Three objects with (J - Kg) < 0.4 are located away 
from the cluster and may be foreground stars. Interestingly, 
LS III -1-46 12 has the second lowest (J - Kg) of all possible 
members, with 0.46 + 0.03, while the location of LS III -1-46 11 
in Fig. 12 shows that it is more reddened (as we already knew 
from the CHORIZOS analysis). We assume that the bulk of clus¬ 
ter members is given by the vertical strip extending between 
(J - Kg) - 0.45 and (J - Kg) - 0.9. This is confirmed by 
their concentration in the IPHAS (r' - i')/{r' - Ha) diagram 



Fig. 6. 2MASS CMD for stars within 3' of the nominal center 
of Berkeley 90 passing the cuts described in the text. The blue 
circles are all the stars with 0.45 < (J - Kg) < 0.9. Red-star 
symbols are used to mark stars with Ha emission, as selected 
from Fig. |7l The thick line is a 2 Ma isochrone with high rotation 
from lEkst^m et al'](l2012l) displaced to DM - 12.0 (log^f = 3.4) 
with E{J - Kg) - 0.75. The brightest objects are LS III -1-46 11 
(to the right of the isochrone) and LS III -1-46 12 (to the left). 

(Fig.llJ, where all except three are distributed in a narrow strip 
with 0.7 ^ (r' - i') % 1.05 that follows the reddening vector, 
and the vast majority have 0.8 ^ (r' - /') 1.0. The main con¬ 

centration of cluster members, lying between LS III h-46 11 and 
LS III h-46 12, shows only a small spread in color, while the stars 
lying immediately adjacent to LS III h-46 11 both to the North 
and West display higher values. The average (J-Kg) for all stars 
with values between 0.45 and 0.9 is 0.67, with a standard devi¬ 
ation cr = 0.12, showing that the objects are evenly distributed 
between these values and therefore defining a typical reddening 
for the cluster is meaningless. 


9 





























J. Mafz Apellaniz et al.: The little-studied cluster Berkeley 90 



Fig. 7. IPHAS (r' - i')l{r' - Ha) diagram. The solid (red) line 
is the locus of the main sequence for £(4405 - 5495) = 1 from 
iDrew et ^ (|2005|) . while the (magenta) dotted line is the locus of 
the main sequence for £(4405 - 5495) = 2, roughly defining the 
envelope for the hot stars that are cluster members. Blue circles 
are used to indicate the stars selected from Fig. and marked 
in the same way there. Since all the stars have been selected to 
have a pseudocolor Qir typical of early-type stars or emission¬ 
line stars and they all have relatively red (J - Ks) colors, we 
can safely assume that stars clearly above the main-sequence lo¬ 
cus are emission-line stars. Their positions in the 2MASS CMD, 
(marked with red-star symbols both here and in Fig. |3, are 
in all cases compatible with this interpretation. Note that both 
LS 111 -1-46 11 and LS 111 -1-46 12 are saturated in IPHAS, as all 
bright O stars in the northern sky are. 


A significant number of sources occupy positions in the 
(r' - i')/(r' - Ha) diagram compatible wit h emission-line stars , 
displaying (r'-Ha) > 0.8 and (r'-T) > 1.2 (ICorradi et al.l2008h . 
Almost all these objects are quite faint and display high values 
of (J - Ks) > 1.4. More than half have Qi^ < -0.1, typical of 
emission line stars. These objects can represent a population of 
pre-main sequence stars associated to the cluster. Interestingly, 
none of them is located in the central region of the cluster. 

As there are no obviously evolved stars, it is not possible 
to assign an age to Berkeley 90. Note that the luminosity differ¬ 
ence between dwarfs and supergiants is relatively small at the 
earliest spectral types. Assuming that LS III -1-46 12 is slightly 
evolved, we can perhaps give an age ~ 2 Ma, but luminous O- 
type stars in Cygnus OB2 with spectral types 03-05 are signifi¬ 
cantly hotter than the main-sequence turn-off dNegueruela et alJ 
|2008|) . By analogy, an age up to 3 Ma is possibl e. For illustra¬ 
tion, F ig. shows a high-rotation isochrone from lEkstrom et akl 
(1201 2h corresponding to an age of 2 Ma, reddened by a repre¬ 
sentative E(J — Ks) — 0.75 and displaced to DM = 12.0 mag 
{\ogd - 3.4). The fit to the position of objects in the central 
concentration is rather good, considering the variable reddening. 
Three objects lie well to the left of the isochrone and therefore 
are either non-members or outliers with low reddening. All of 
them lie to the East of LS III -1-46 12, except one that is located 
at the northern edge of the region analyzed. However, removing 
these objects does not change the average color or its standard 
deviation. On the other hand, a few objects lying to the right of 
the isochrone, with (J - Ks) ~ 1.0 are likely cluster members 
with higher than average reddening. 


In spite of the presence of three early O-type stars (two in 
LS III -1-46 11 and at least one in LS III -1-46 12), Berkeley 90 
seems to contain very few OB stars. Only one photomet¬ 
ric member is sufficiently bright in Ks to be a late-O star, 
2MASS J20350798-H4649321 (Eig.[B, and this object has a po¬ 
sition in the IPHAS diagrams consistent with being an emission¬ 
line star (it is the object marked with both a blue circle and a 
red-star symbol in Eigs. |6] and |7]i. The bulk of the population 
starts almost 3 magnitudes below LS III -1-46 12, at Ks ~ 10.5, 
an intrinsic magnitude roughly corresponding to a B1 V spectral 
type. 


4. Discussion 


O stars earlier than type 04 are very scarce in the Galaxy. 
Prior to this work, there were only two examples known in 
the northern hemisphere, Cyg OB2-7 and Cyg OB2-22 A 
(IWalborn et al.112003 ISota et al1l201 ih . Given that we know so 
few very massive stars, it is crucial to keep searching for them to 
increase our statistics if we want to establish what is the stellar 
upper mass limit and its dependence on metallicity and environ¬ 
mental conditions. 

We have shown that LS III -1-46 11 as a very massive eccen¬ 
tric binary composed of two near-twin stars. A few years ago 
this may have been seen as a fluke but several similar systems 
have been discovered recently. 


HD 93 129 Aa Ab _ (iNelan et all l2^ I 2 OIOI: 

iMaiz Anellaniz et ^ 120051 2008; Sota et alF^ldl) may 
be even more massive, given that the primary is of spectral 
type 02 If*, but the secondary is one magnitude fainter and 
the orbit is yet undetermined and appears to be decades to 
centuries long; __ 

Cyg OB2-9 (INaze et aUlTOTl 05-5.5 I -H 03-4 III) has a 
very similar mass ratio and a somewhat larger eccentricity, 
its (Ml -I- M 2 ) sin^ i is only slightly lower, its period is an 
order of magnitude larger, and the spectral types are slightly 
later. __ 

R139 (iTavlor et alJl201 ll 06.5 lafc - 1 - 06 laf) is also similar 
to LS III -1-46 11 in terms of period, eccentricity, and mass 
ratio but the two stars are mid-O supergiants and their lower 
mass limits are significantly higher, between 60 and 80 Mq. 
Those masses indicate that not all objects above 60 Mq are 
WNh stars (of course, as long as we do not know the inclina¬ 
tions in this and other cases, we will not know what the true 
masses are). 

Two additional examples of very massive stars in 
elliptical orbits whose large eccentricity made the 
discovery of t heir binarity diffic ul t are WR 22 s 
HD 92740 


119791 LSchweickhardtet al.1 119991 WN7 08-9.5 III:) 

and HD 93162 (Gamenetal. l2008t LSota et all 12014 

02.5 If*/WN6 - 1 - OB). 
There are also very mass 
circular orbits such a 
l2004 ISchnurr et al.l 

ive twin systems in shorter, near- 
s NGC 3603-A1 (iMoffat et al.l 
20081 WN6ha - 1 - WN6ha), 

WR 20a (iRauw et al 

2004 iBonanos et al.1 l2004 

ICrowther & Walborn 20 

3 03 If*AVN6 -1- 03 If*AVN6), 


and Pismis 24-1 (iMafz Anellaniz et al.1 120071 
03.5 If* -I- 04 Ill(f) -I- ...). The latter also includes a 
third very massive component in a long orbit. 

All of those systems are located in regions with large num¬ 
bers of O stars (Carina Nebula; Cygnus OB2; 30 Doradus; 
NGC 3603; Westerlund 2; and, to a lesser extent, Pismis 24). 
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LS III -1-46 11 is the oddity in that respect, being located 
in a signihcantly less massive cluster or association. In 
that respect, a more similar case may be HD 150136, 
which is apparently less evolved but is a triple system 
[03-3.5 V((f*)) -I- 05.5-6 V((f)) -h 06.5-7 V((f))] with 
a most massive star of 53+10 M© in a relatively small 
cluster, NGC 6193, with another nearby star, HD 150135 
[06.5 V((D) zl, yielding a similar makeup to that of 
Berkeley 90 (iNiem ela & Gamed 1 ^051: Mahv et al. 20121: 


^na et alJl2013btlSanchez-Bermudez et al.ll2013 : Sota et all 

2014 . 


All of the cases mentioned above refer to distant very mas¬ 
sive stars (Cyg OB2-9 and HD 150 136 are the closest ones but 
they are beyond 1 kpc) but what should be even more surprising 
is that some of the brightest and closest O stars in the sky have 
been recently discovered to have e ccentric companions. That is 
what has happened with 0' Ori C ( Krau^_e^^ 2007^ Sota et all 
1201 iL 07 Vp + ...) and cr Ori A dSimon-Dlaz et alJ2()l iL 2015aL 
09.5 V + BO/1 Vn, note that sigma Ori B is B0.5 V). Why 
have not these systems been discovered before? There are two 
reasons: eccentric systems require extensive spectroscopic mon¬ 
itoring since their velocity differences may be too small to be 
detected during a large fraction of their orbits (as it happened 
with LS III +46 11) and in some cases interferometry is the only 
way to detect their multiplicity because of their large semi-major 
axes. It has not been until the last decade that large-scale spectro¬ 
scopic monitoring of many O stars has started and that interfero¬ 
metric technology has allowed similar surveys using those tech¬ 
niques. However, many systems still remain outside the reach 
of such surveys so discoveries should continue in the following 
years. 

Our data cannot provide conclusive results on the masses 
of LS III +46 11 and LS III +46 12. Without eclipses, we can¬ 
not accurately measure the inclination of the LS III +46 11 orbit 
and ourkeplerian masses of 38.80+0.83 M© and 35.60+0.77M© 
have the sin^ i factor included in them, making them just lower 
limits. Note that the star is not prese nt in the public relea se of 
the Northern Sky Variability Survey (IWozniak et al.ll2 004 and 
that the time coverage in the Super WASP ( Pollacco et al.ll20()^ 
public archive is very limited, so a thorough search for eclipses 
is not possible at this time. In any case, eclipses are unlikely, 
since the large separations in the orbit would require an incli¬ 
nation very close to 90 degrees in order for them to take place. 
The inclination could be constrained in the future with broad¬ 
band polarimetry or through the phase-dependent behaviour of 
excess emission from colliding winds 

For evolutionary masses, we have a different problem: our re¬ 
sults for the spectral types and distances are inconsistent. There 
are three possible explanations: 


1. A straightforward interpretation of the CHORIZOS results 
place LS III +46 11 at a log^f between 3.45 and 3.67 (in 
parsecs, after correcting for the existence of two stars) and 
LS III +46 12 at a logd between 3.10 and 3.28. Therefore, 
one possibility would be that both objects are physically un¬ 
related and just the result of a chance alignment. In this case, 
the luminosity classihcation criteria for early-type O stars 
would retain their physical meaning and there would be no 
need to invoke the existence of undetected multiple systems. 
Nevertheless, we judge this possibility to be highly unlikely, 
given the small number of early-type O stars that exist and 
the existence of the underlying cluster Berkeley 90. In any 
case, this hypothesis will be tested soon once the Gaia paral¬ 
laxes become available. 


2. An alternative would place both LS III +46 11 and 
LS III +46 12 at a logt/~3.4, consistent with the 2MASS 
photometry of Berkeley 90, but would require that the cur¬ 
rent luminosity classihcation criteria for O stars based on 
the depth of Hen 44685.71 does not really rehect a func¬ 
tion or luminosity but instead is just a measurement of wind 
strengtlQ. We consider this option unlikely, as there are 
both theoretical reasons why (for the same Teff and metal- 
licity) wind strength should strongly correlate with lumi¬ 
nosity and observational data t hat corroborate that associ¬ 
ation (e.g. IWalborn et af]l20l4 . This hypothesis could be 
tested by obtaining a good-S/N high-resolution spectrum 
of LS III +46 12 and modelling it with e.g. FASTWIND 
or CMFGEN. Those codes derive gravity from a ht to the 
Balmer line prohles which is independent of He ii /14685.7L 

3. A third option is that LS III +46 12 is another near-twin 
binary system composed of two very-early-type O dwarfs. 
In this scenario, logd would be in the 3.40-3.45 range, 
marginally consistent with the spectroscopic parallaxes for 
LS III +46 11 and Berkeley 90. This solution may be ad hoc 
but we believe it to be the most likely one. Indeed, there is 
a precedent with an object that is a near-s pectroscop ic twin 
of LS III +46 12, HD 93 250 (04 III(fc). lSotaetaHlM4 . 
Its spectroscopic parallax was incompatible with the well 
known value of the Carina Nebula u ntil it was discov ered to 
be a binary through interferometry ([Sana et alJl201 Ih . Note 
that there is a large range of separations for which we would 
not detect signihcant velocity variations if the system were 
a spectroscopic binary, that the inclination could be small, 
that our spectroscopic campaign has not been as thorough 
for LS III +46 12 as it has been for LS III +46 11, and that 
X-ray excesses due to wind-wind collisions are expected to 
be weaker for dwarfs than for supergiants. 

Related to the issue of the masses is the age of Berkeley 90. 
Until one of the above explanations is conhrmed, we cannot 
provide a final answer. We should note, however, that if the 
third option is the correct one, it is possible to derive a con¬ 
sistent age of 1.5-2.0 Ma for LS III +46 11, LS III +46 12, and 
Be rkeley 90 using the low-rot ation Geneva evolutionary tracks 
of iLeieune & Schaer^ (1200 Ih . Under those assumptions, the 
masses of the two stars in LS III +46 11 would be in the 60- 
70 M© range and those in LS III +46 12 would be 35-45 M©. 

The results in this paper reveal that Berkeley 90 can be an 
important cluster to resolve the issue of how the stellar Initial 
Mass Function (IMF) is sampled. There are two alternative the¬ 
ories, one that states that the IMF is built in a sorted way, with 
many low-mass stars being formed before any high-mass star 
can appear, and another one that states that the IMF is built 
in a stochastic way, with some rare cases where it is possible 
to form high-m ass stars with a relat ively small number of low- 
mass ones (see iBressert et al.ll2012l and references therein). If 
the third option on the LS III +46 11 and LS III +46 12 masses 
above is correct, we would have a cluster with four stars above 
35 M©, one of them with ~70 M©, and none or just one 
(2MASS J20350798+4649321) between 20 and 35 M© (see 
Fig. |6] and the associated discussion in the text). Using the total 
stella r mass above 20 M© and assuming a Kroupa IMF (iKrounal 
120011) between 0.1 M© and 100 M©, we derive that the mass 


’’ Also note that the continuum of LS III +46 12 to the right of 
Hen 44685.71 is slightly rised, possibly signalling an anomalous pro¬ 
file indicative of luminosity class III rather than V, which would have a 
similar effect to the one described here. 
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of Berkeley 90 is =s2000 M©, which is consistent with the ap¬ 
pearance of the cluster in relationship with other similar objects. 
However, such a large maximum stellar mass in such a small 
cluster could be incompati ble with the sorted sampling scenario 
dWeidner & Kroup3l2006h but can be accomodated within the 
stochastic scenario: a 2000 M© with the conditions above should 
have, on average, 2.3 stars in the 20-35 M© range and 1.6 stars 
in the 35-100 M© range. We observe 1 and 4, respectively i.e. 
off from the expected values but within the range of possibili¬ 
ties expected for a single case within the number of similar clus¬ 
ters in the solar neighborhood. We need better data (constraining 
the nature of LS III -1-46 12 with interferometry, studying a deep 
CMD of Berkeley 90 to accurately measure the cluster mass, and 
resolving the distance issues with Gaia) to provide a hnal answer 
but Berkeley 90 appears to be a good candidate for a small clus¬ 
ter with a star too massive for the sorted sampling scenario. 


5. Conclusions 

- Berkeley 90 is a young stellar cluster dominated by two early 
O-type systems. LS III -1-46 11 is an SB2 composed of two 
very similar 03.5 If* stars. LS III h- 46 12 is spectroscopi¬ 
cally single and has a spectral type 04.5 V((f)). 

- LS III -1-46 11 has an eccentric orbit (e between 0.56 and 
0.57) with a 97.2-day period and minimum masses of 
38.80+0.83 M© and 35.60+0.77 M©. Since we do not 
know the inclination we cannot calculate accurate keplerian 
masses. 

- LS III +46 11 has a significantly higher extinction than 
LS III +46 12. The optical+NIR extinction law is close to 
the average one in the Galaxy. 

- There are no apparent bright visual companions to either sys¬ 
tem. 

- The evolutionary masses of LS III +46 11 andLS III +46 12 
are incompatible with the two systems being located at the 
same distance and having the same age. We consider dif¬ 
ferent solutions to the problem and consider that the most 
likely one is the existence of an undetected companion for 
LS III +46 12, for which there is plenty of room in terms of 
period, inclination, and eccentricity not yet explored. 

- Berkeley 90 is a cluster with considerable differential extinc¬ 
tion and its stellar mass is possibly too low to harbor both 
LS III +46 11 and LS III +46 12 under the sorted sampling 
scenario for the IMF. 
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